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Abstract 

Slow and accurate finger and limb movements are essential to daily activities, but their neural control and 
governing mechanics are relatively unexplored. We consider neuromechanical systems where slow movements 
are produced by neural commands that modulate muscle stiffness. This formulation based on strain-energy 
equilibria is in agreement with prior work on neural control of muscle and limb impedance. Slow limb 
movements are driftless in the sense that movement stops when neural commands stop. We demonstrate, in 
the context of two planar tendon-driven systems representing a finger and a leg, that the control of muscle 
stiffness suffices to produce stable and accurate limb postures and quasi-static (slow) transitions among 
them. We prove, however, that stable postures are achievable only when muscles are pre-tensioned, as is the 
case for natural muscle tone. Our results further indicate, in accordance with experimental findings, that 
slow movements are non-smooth. The non-smoothness arises because the precision with which individual 
muscle stiffnesses need to be controlled changes substantially throughout the limb’s motion. These results 
underscore the fundamental roles of muscle tone and accurate neural control of muscle stiffness in producing 
stable limb postures and slow movements. 
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Introduction 

Fine grasp and manipulation of objects such as threading a needle or writing require accurate finger 
and limb postures and slow movements. Inaccuracies in performing such tasks often arise in neuromuscular 
conditions such as spasticity, Parkinson’s disease and stroke HIES]. An understanding of these inaccuracies 
and their causes necessitates a deeper knowledge of the neuromechanical conditions required for the muscles 
to produce slow and accurate limb movement. The neural control of muscles is complex [22]; however, it 
is believed that the neural input determines the force and stiffness of the muscle [18, [19; [36j [37l 41, EL. 
We therefore propose a computational framework for slow movements where the neural control of muscle 
activity is simulated by altering the muscle stiffness properties. 

Slow limb movements belong to a class of mechanical systems, known as driftless , where movement stops 
when actuation stops. Drift less systems arise in many applications, including satellite dynamics mug, 
robotic vehicles [14 , and biolocomotion Ha Eg Eg Eg. We apply these principles to examine the neurome¬ 
chanical properties of quasi-static limb movements. This viewpoint is reminiscent of but not identical to 
the equilibrium point hypothesis (EPH), which states that limb movements are generated by a sequence of 
equilibrium points along a desired trajectory mmm- EPH considers that the magnitude of the “control” 
force exerted on the limb, at any time, depends on the difference between the “actual” limb dynamics and 
the desired equilibrium point, and derives a second-order differential equation for the “error” between the 
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actual and desired dynamics mmm- Here, we are not concerned with the relaxation dynamics. We posit 
that this relaxation time scale is much faster than the time scale associated with slow limb movements, and 
the limb reaches a state of equilibrium instantaneously in response to changes in muscle stiffness. We seek 
to determine whether equilibrium postures of multi-joint limbs and quasi-static transitions among them can 
be achieved by modulating the muscles’ stiffness values. 

We demonstrate that the control of muscle stiffness suffices to produce stable limb postures and quasi¬ 
static transitions among them. This is in agreement with results from impedance control [23] GS]* However, 
we go beyond these results to prove that stable postures are achievable only when muscles are pre-tensioned 
as in the case of natural muscle tone. Muscle tone is defined clinically as a state of active muscle contraction 
at rest [43] . and depends on the physiological or pathological state of the limb. Healthy muscles at rest are 
always producing a baseline level of force at the tendon [58]. This baseline force is equivalent to the muscle 
“threshold” value of EPH [17]. Unlike threshold control, we do not alter the pre-tensioning length directly. 
Rather, we define muscle tone as a baseline non-zero tension in the tendon, and view muscle stiffness as a 
neurally controlled parameter. 

Controlling muscle stiffness to produce equilibrium limb postures does not admit a unique solution 
due to the over-actuated nature of the tendon-driven systems, thus, the same posture can be achieved at 
multiple combinations of stiffness parameter values. Each combination has a distinct strain energy cost. 
If a relatively-low energy level is inaccessible, say due to a pathological condition, the same posture may 
be realized at a higher energy cost. Therefore, we formulate the problem of muscle stiffness control as an 
optimization problem [25, 39] that minimizes the strain energy of the limb. Our optimal stiffness values 
are in agreement with the reciprocal-inhibition phenomenon: relaxation of the antagonist muscle during 
activity of the agonist [46] . Further, we find that, for constant joint moment arms and muscle pre-tensioning 
values, optimal stiffness values depend on the change in the joint angle but not on the reference posture. 
This symmetry with respect to reference postures is broken in the case of the anatomically more realistic, 
posture-dependent moment arms. 

Finally, we investigate the smoothness of quasi-static transitions among limb postures given physiological 
limits on the neural control of muscle stiffness. Discontinuities, or non-smoothness, in slow limb movements 
have been reported in experimental studies [S3IM1ES], and are thought to be due to non-trivial coupling 
among multiple neural and peripheral sources mmmm- The intrinsic stiffness of individual muscles 
is linked to muscle activation and force [12] [24] [9:. Therefore, the precision with which muscle stiffness 
can be controlled is in the same general range as the precision with which muscle force can be controlled. 
Prior work [201 EH EZl 121] shows that, at low force productions, muscle force can be regulated with a 
precision in the range of 2.5 to 10%. Under such limitations, we observe discontinuities in the limb’s quasi¬ 
static trajectories, indicating unreachable postures because of insufficient precision in the muscle stiffness 
parameters. Our results indicate that these discontinuities are spatially localized in the limb workspace. 


Methods 


We examine the slow movements of tendon-driven limbs in the context of two planar model systems: an 
idealized finger model and a more anatomically-realistic leg model. For notational convenience, we consider 
the general case of a planar limb of n joints and let /i,..., l n be the lengths of the individual limb segments. 
The limb posture is defined by the column vector of joint angles 6 = (0 1 ,..., 0 n ) T where the superscript () T 
denotes the transpose operator. The joints are driven by m muscles (m > n ) of variable stiffness parameters 
ki , i = 1,..., m, whose tendon paths produce an n x m moment arm matrix R(0) [49 , 
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(1) 


Each entry (here, i = 1,..., n, j = 1,..., m) denotes the “rotation” the j th tendon produces at the i th 
joint. It is positive when pulling the j th tendon induces a positive rotation (i.e., counterclockwise rotation 
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Figure 1: Finger (a,b) and leg (c,d) models: kinematics (a,c) and muscle routing (b,d) overlaid on top of the physiological 
system. 


per the right-hand rule). Depictions of the finger (n = 2) and leg (n = 3) models and their tendon paths are 
shown in Figure [l] Note that the extensor mechanism of the fingers is rather complex [SB, and Figure [ljb) 
corresponds to an idealized case of multi-articular muscles whose moment arms change signs. The tendon 
paths of the leg model in Figure [ljd) are more anatomically-realistic [SI ESI E]- For clarity, Figure m 
shows only seven of the fourteen muscles/muscle groups, amounting to a total of twenty-three individual 
muscles that we use in this leg model. In both the finger and leg models, we examine two types of joints: 
simple hinge joints where the moment arm matrix R is constant for all values of 0, and non-circular joints 
where the moment arm matrix is posture-dependent. In the finger model, the constant moment arm matrix 
is chosen to be 

/ -0.8 -0.8 0.7 0.7 \ 

K_ v “°- 5 °- 2 “°* 5 °- 2 ) ’ 

in length units of cm. These values lie in the range of average moment arm values for the index finger 
reported in [U|52]. For the posture-dependent moment arms, we define joint ellipses with semi-major and 
semi-minor axes in the range of measured index moment arms found in pQ. In the leg model, we utilize the 
posture-dependent moment arm values of the model developed by [2 . 

A change in limb posture corresponds to a rotation A0 = 0i — 0o of the joints from a reference limb 
posture #o to a new limb posture 6\. A rotation A0 fully determines the excursions As = (Asi,..., As m ) T 
of all muscles [DEO], 

pO i 

As = - / (R(0 )) t dd = -R t A e , (2) 

Je o 

where the negative sign indicates that a positive rotation of the joint will shorten the muscles that in¬ 
duce it; and vice versa. For constant moment arms, the second equality in © is evident. In the case of 
posture-dependent moment arms, because we only consider small posture changes A0, we use the trapezoid 
integration rule such that equation ([2| holds with R = (R(0 q) + R(0 1 ))/ 2 , where R(0o) and R(0i) are the 
moment arm matrices evaluated at postures Oq and 6 1 , respectively. 

We assume that each muscle acts as an ideal spring whose stiffness parameter is set by the neural 
command isuiTi mus] . We also allow each muscle to be pre-stretched from its resting length by an amount 
Al 0 = (AZi,..., A lm) T . The total strain energy of the system, E , in the presence of muscle excursions As 
is then given by 


E = I (As + A1 0 ) t K(As + Al 0 ) . 


(3) 
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As per the impedance, in this case stiffness, control formulation [22], the stiffness matrix K is an m x m 
diagonal matrix of entries ki. 


K = 


fki 

Vo 



(4) 


The corresponding forces f = (/i,..., / m ) T generated by the muscles are given by 


f = -K(As +Al 0 ) • (5) 

These muscle forces, in turn, produce torques r = (ti, _, r n ) T at the joints defined by r = Rf. A stable 

posture is achieved by satisfying the static equilibrium condition at the joints, namely, the total torque at 
each joint must be zero, r = Rf = 0. Slow, i.e., quasi-static, limb movements are achieved by sequentially 
satisfying this static equilibrium at each posture. 

In the absence of external loading, if the muscles are not pre-tensioned, setting Al 0 = 0 and substitut¬ 
ing © and ([5| into Rf = 0 gives 

RKR T A0 = 0 , (6) 

which only admits the trivial solution AO = 0 for all K because the n x n square matrix RKR T is full rank 
when R is full rank. Thus, there exists only one solution at the reference posture, and other stable postures 
are not achievable by varying K. In contrast, when the muscles are pre-tensioned, one has 

-RK(-R T A0 + Alo) = 0 . (7) 


The system is thus unconditionally controllable and equilibrium postures are achievable by proper choice of 
muscle stiffness K. 

The quasi-static formulation implies that, starting from a given initial posture, the limb has to transition 
to “nearby” equilibrium postures, that is to say, A0 must be small. A large value of A0 means that the limb 
would instantaneously jump from its current posture to a far away equilibrium posture, which violates the 
slow movement assumption. Therefore, in all subsequent analyses, we consider quasi-static transitions to 
only nearby postures. In particular, we formulate: (I) an optimality problem where for a given equilibrium 
posture 0 , we solve for optimal stiffness values that minimize the strain energy function associated with 
transitions to nearby postures, and (II) a reachability problem where we explore reachable equilibrium 
postures for given stiffness values. 


Optimal equilibrium postures. Starting at a posture 0o, a desired limb posture 0 near 0 q can be achieved 
by tuning the muscle stiffness parameters so that the vector (—R T A0 + A1 q) in equation ([7]) lies in the 
null space of the n x m matrix RK. The resulting stiffness values are not unique, in general. For each 
limb posture 0, equation ([7]) yields a family of solutions K, at multiple strain energy levels, for which 0 
is achievable. This redundancy leads us to look for optimal stiffness values K opt that minimize the strain 
energy function while satisfying the equilibrium constraint and inequality constraints on permissible muscle 
stiffness values, namely, 


min 

K 


E = i(-R T A 9 + A1 o ) t K(-R t A 0 + Al 0 ) 


subject to — RK (—R T A0 + Al 0 ) = 0 , 


and 


<ki<k 


max • 


(8) 


Reachable equilibrium postures. Given a specific combination of muscle stiffness parameters K, the achiev¬ 
able equilibrium posture 0 can be found by solving the “forward problem” in ([7]), which we rewrite as 
RKR t A0 = RKAIq. Given that the matrix RKR T is invertible, solutions to the forward problem are 
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unique: for each K, there exists one and only one A0, which in turn gives a uniquely-defined new posture 
0. As in the optimality problem, we only examine small changes of equilibrium posture AO. 

A few remarks on the parameter values of A: m ^ n and kmax are in order here. k m in is always greater than 
zero ensuring that the optimization problem is well-posed. Whereas leg muscles generate larger forces than 
finger muscles eh m ns], the muscle stiffness values are comparable in both, with typical values of a few 
hundred Newtons per meter (N/m) (42 ] |T0l 13 SO] • Without loss of generality, we choose a range of muscle 
stiffness values kmin = 100 N/m to k max = 1000 N/m, in agreement with these experimental data. We 
now non-dimensionalize all parameters using a characteristic length scale Z* = 10 cm and a characteristic 
force /* = 100 N. Non-dimensional parameters enable us to obtain generic results, which we can later 
easily scale up or scale down to different limb sizes in both generic and subject-specific modeling. Using the 
characteristic length and force scales, the dimensionless range of stiffness values is k m i n = 0.1 to k max = 1. 
The normalized values of the pre-stretched lengths A^ are all set to 0.1 in the finger model whereas to 0.5 
in the leg model. Note that in addition to constant pre-stretch, we consider pre-stretched lengths that vary 
with posture, as described in the Results section. 


Results 

Optimal equilibrium postures. We consider first the simpler finger model since it is more amenable to visual 
representations. Given 0 m (0i, 02) T , the finger’s posture in the (x,y )~plane follows in a straightforward 
way. Its workspace W is the set of all points ( x w ,y w ) of the (x,y )~plane accessible by the finger’s endpoint 
for all admissible 0i and # 2 , 

x w = h cos#i + l 2 cos(0 2 — 0 i), y w = h sin 6b + l 2 sin(0 2 — 6b). (9) 

A depiction of the joint angles space and the workspace of the finger model is shown in Figure |2j Here, the 
admissible joint angles are taken to lie in the range 6b E [—80°, 10°], 0 2 E [0°, 90°]. Initial postures Oo in this 
admissible range are discretized using a 9 x 9 regular grid, marked by x. For each 0 0 , we explore optimal 
transitions to equilibrium postures in a small square neighborhood centered at 0 0 and of side length 10°. 
More specifically, we consider a total of nine nearby equilibrium postures, including staying at the initial 
posture Oo. We compute the optimal stiffness values for the change in joint angles AO associated with a 
transition from the initial posture to each of these nine postures. We linearly interpolate these values in the 
neighborhood of each initial posture to construct a complete map of optimal stiffness values. 

Figure [3] depicts the optimal stiffness values K opt for the finger model with constant moment arms. 
Muscles m 2 and m 4 exhibit higher stiffness values than m 1 or m 3 , which means, considering that stiffness 
values are proxy to the intensity of the neural command to the muscle, that muscles m 2 and 777,4 require 
greater neural command. A closer look at the stiffness variation within each neighborhood shows that, in 
the context of constant moment arms, all neighborhoods are identical, that is to say, optimal stiffness values 
depend only on the change in joint angles A0, or posture change, and not on the initial posture. Further, 
the stiffness variation within each neighborhood is such that m 2 and 777,4 counter-balance each other in the 
0i-direction - as the stiffness of m 2 increases, that of 777,4 decreases and vice-versa. This is consistent with 
the notion of agonist-antagonist muscle pairs. If one muscle acts as an agonist and contracts to rotate the 
joint, the torque that it produces must be counterbalanced by another muscle to prevent the system from 
accelerating. While this language of co-contraction of agonist-antagonist pairs does not extrapolate well to 
multi-muscle, multi-joint systems m, the ideas that muscles produce postural equilibrium via co-contraction 
are well established mM- 

Figure [4] depicts the optimal stiffness values for the case when the moment arms vary with joint angles. 
As in the case of constant moment arms, muscles m 2 and 777,4 require higher stiffness values, and the 
neighborhood around each initial posture indicates that each muscle has a “preferred” direction of higher 
stiffness, which is the same as in Figure [3J However, the optimal stiffness values now depend on moment 
arms, and thus on initial posture. Interestingly, the stiffness variation across the whole workspace is also 
characterized by the muscles m 2 and 777,4 counter-balancing each other, indicating a hierarchical structure 
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Figure 2: Finger model: (a) configuration space and (b) the corresponding endpoint space which deforms the squares into 
” diamonds.” 



Figure 3: Finger model with constant moment arms: optimal stiffness parameters needed to achieve each posture in (a) 
configuration (joint angle) space, and (b) endpoint space. The plots show the activation level of each muscle mj, i = 1,..., 4, 
that minimize the strain energy function. Note the agonist-antagonist actions of m 2 and 1714 in the 6\ -direction. 


in the way this agonist-antagonist muscle pair acts in transitions to nearby equilibrium postures and to 
postures across the whole workspace. 

In the leg model system, initial postures are taken in the range 0\ G [—20°, 100°], O 2 G [0°,100°] and 
0s G [—20°, 30°], and twenty-seven equilibrium postures are computed in the vicinity of each initial posture. 
Figure |5|a) shows the optimal stiffness combinations in joint angle space for constant moment arms, while 
Figure pFb) for variable, posture-dependent moment arms. For clarity, only the leg muscles that manifest the 
highest stiffness values are shown, namely, muscles gluteus maximus , rectus femoris , tensor fasciae latae and 
tibialis anterior , labeled mg, 77717 , 77721 and 77122 , respectively. The remaining muscles, while they individually 
have low stiffness values, contribute collectively to the equilibrium postures. As in the case of the finger 
model, optimal stiffness combination K opt depends on the change AO but is independent of the initial joint 
angles in the case of constant moment arms but depends on both in the case of posture-dependent moment 
arms. 
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Figure 4 : Finger model with posture-dependent moment arms. Note the agonist-antagonist actions of m2 and 7714 in the 
61 -direction, both locally and over the whole space of joint angles. 




Figure 5 : Leg model: optimal stiffness parameters needed to achieve each posture in (a) configuration (joint angle) space,and 
(b) endpoint space. The plots show the activation level of the four highly active muscles, ms 77117, 77121 and 77122 (gluteus 
maximus , rectus femoris , tensor fasciae latae and tibialis anterior , respectively), that minimize the strain energy function. 


Quasi-static trajectories. We consider quasi-static trajectories where the limb endpoint is required to slowly 
trace, back and forth, curved and straight lines in the workspace as depicted in Figure [6] We solve for 
the optimal muscle stiffness combinations required by the limb to perform this back and forth motion. To 
this end, the quasi-static trajectories are discretized, and the optimization problem in equation is solved 
sequentially along the discrete trajectory. 

We investigate the effect of constant versus variable muscle pre-tensioning. Constant pre-tensioning 
means that the system is “memoryless” because muscle pre-tensioning does not depend on the previous 
muscle configuration. Variable pre-tensioning is defined according to the update rule depicted in Figure [7| 
where the pre-tensioning value at each posture is equal to that at the previous posture plus the change 
in muscle length that occurred as a result of moving from the previous to the current posture; thus the 
pre-tensioning has “memory.” 
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Figure 6 : Endpoint trajectories for (a) and (b) finger model and (c) leg model. 
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Figure 7 : Variable pre-tensioning update rule: the pre-tensioning is history-dependent, taking into account the muscle length¬ 
ening (a) or shortening (b) that took place in the previous posture change. 


Figure [8] depicts the optimal stiffness values for the finger model tracing the trajectories of Figure j6^a,b) 
for all combinations of moment arms and pre-tensioning scenarios, whereas Figure [9] shows similar results for 
the leg model tracing the curve of Figure |6jc). In both models, the optimal stiffness values are characterized 
by a local jump as the limb reverses its motion to trace the trajectory backward. The jump in the stiffness 
of certain muscles gets attenuated in the case of variable moment arms, but all jumps get smoothed out 
completely when the muscle pre-tensioning is changing continuously with posture. In these cases, one 
observes the relaxation of the antagonist muscles during the activity of the agonist, and the opposite effect 
is clearly seen during the return portion of the movement. For instance, the gluteus maximus (mg) and 
rectus femoris ( 77117 ) in the leg model act against each other to complete the trajectory. 

Reachable equilibrium postures. We explore the reachability or feasibility of tracing certain trajectories 
within the limb’s workspace given the inherent constraints on the physiological precision in the neural 
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Figure 8: Finger model: optimal stiffness values for the quasi-static trajectories shown in Figure [6ja,b). 
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Figure 9: Leg model: optimal stiffness values for quasi-static curved trajectory shown in Figure [SJc). The most active muscles 
during this trajectory are the gluteus maximus (rnj and ms), rectus femoris (m 17 ), biceps femoris long head (mig), tensor 
fasciae late (m 2 i) and tibialis anterior ( 77122 )- 


control of muscle stiffness. The precision with which muscle activation can be controlled is not infinite; 
it is limited by the mechanisms of recruitment and rate coding used to control muscles. Thus, the space 
of possible stiffness values K is not infinitely smooth. Here, we systematically sample the stiffness space 
at different precision levels, i.e., degrees of discretization granularity Afc. For each precision level Afc, we 
calculate the set of admissible departures A0 from a reference posture, and consequently the set of reachable 
(feasible) postures 0. Our goal is not to minimize the strain energy at each reachable posture, but to find 
the set of nearby reachable postures as a function of the precision with which individual muscle stiffnesses 
can be controlled. Minimizing the strain energy at each reachable posture introduces additional constraints, 
making it even more challenging to reach a desired posture at a given stiffness precision level. 

We discretize the range of stiffness values from fc min to k max using constant increments Afc, producing 
M = (kmax — kmin )/Afc + 1 discrete stiffness values for each muscle. For m muscles, this discretization 
generates M m different stiffness matrices K, which amounts essentially to a uniform sampling of the Tri¬ 
dimensional muscle stiffness space. For example, setting Afc = 0.05 in the finger model amounts to 19 4 
distinct K matrices. For each K, we solve equation ([7]) to compute A 0 starting from a reference posture 0 $ , 
and use © to compute the reachable endpoint location (x w , y w ). We constrain the full set of M m reachable 
points to those that are in the vicinity of the reference posture. 

For illustration purposes, we investigate the ability of the finger model to trace multiple trajectories for a 
given Afc. Starting from the initial posture, we sequentially transition to the next point along the trajectory 
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Figure 10: Finger model: achievable versus desired trajectories overlaid on top of the reachable postures. The level of stiffness 
precision affects the reachability. As precision increases (from left to right), the endpoint trajectories become smoother and 
closer to the desired trajectories. Discontinuities in the set of reachable points vary nonlinearly in the workspace, e.g., regions 
highlighted in red. 


by locating the nearest admissible reachable posture. Figure [To] shows the desired versus reachable endpoint 
trajectories for three precision levels A k = 0.1,0.08,0.05, superimposed on the set of all reachable point 
in the limb workspace. The reachable paths show discontinuities. These discontinuities decrease as the 
resolution of muscle stiffness increases (going from A k = 0.1 and A k = 0.05). The discontinuous behavior 
arises from unreachable desired postures, in which case the limb settles for the nearest admissible posture 
from the set of reachable postures. Further, we observe that discontinuities are spatially localized, implying 
that the demands on the precision in muscle stiffness varies nonlinearly in the limb workspace. 


Discussion 

We developed a novel mathematical formulation for slow and accurate limb movements using the frame¬ 
work of driftless mechanical systems, and considered the neural control of muscles to change the muscles’ 
stiffness parameters. We used this framework to address three fundamental questions in slow limb move¬ 
ments: (1) can slow limb movements be achieved by modulating the muscle stiffness parameters only? (2) 
what is the role of muscle tone? and (3) how are these movements affected by limits on the precision with 
which muscle stiffness can be controlled? We demonstrated that muscle stiffness control is sufficient to 
produce accurate and slow limb movements, but only when the muscles are pre-stretched, as in the case 
of natural muscle tone. We probed the limb’s ability to trace prescribed trajectories under constraints on 
muscle stiffness precision. We identified discontinuities in trajectory tracking. By comparing different reso¬ 
lutions of muscle stiffness, we found that high resolution of muscle stiffness is key for the limb to improve 
its accuracy in tracking a trajectory. 

The discontinuous nature of slow finger movements was reported in numerous experimental studies 
Ei mg. Several studies attributed these discontinuities to neural sources while others held responsible pe¬ 
ripheral stretch reflexes [Sol USUIS]. Recent reports seem to suggest an inter-related network that includes 
both the brain and the periphery naim. Our work supports the latter view of a combined neuro-mechanical 
origin of the discontinuities and presents an alternative and novel explanation of the mechanism mediat¬ 
ing these discontinuities, namely, the precision in muscle stiffness. Muscle stiffness is known to have two 
components: (i) short-range or inherent stiffness of the muscle coming from sarcomere mechanics, motor 
unit recruitment and rate coding, and tissue properties; and (ii) reflex-driven stiffness coming from the 
stretch-reflex response mediated by spinal and transcortical afferent circuits mmu- Pathological dis¬ 
ruption of either mechanism could degrade the level and precision of muscle tone and muscle stiffness. One 
example is abnormally high muscle tone, or hypertonia, which is defined as abnormally increased resistance 
to externally imposed movement about a joint, and may be caused by spasticity, dystonia, rigidity, etc. [43] . 
Other examples include pathologic tremor and bradykinesia (slow movements) in people with Parkinson’s 
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disease mmm, and stroke [33]. Our results suggest that degradation in the precision of muscle stiff¬ 
ness would exacerbate the inherent non-smoothness of slow movements in tendon-driven limbs and provide 
alternative and complementary peripheral explanations for the aforementioned clinical symptoms. 

We made a number of simplifying assumptions. The assumption that small finger muscles and larger 
lower limb muscles have similar stiffness range is based on recent reports that relatively small finger mus¬ 
cles, such as the extensor digitorum or flexor carpi radialis , have stiffness values m similar to a larger leg 
muscle, such as the rectus femoris (around 3007V/ra at rest) [40] . We neglected gravity in both the finger 
and the leg model but we obtained similar results in the presence of gravitational effects (not shown for 
brevity). We considered the initial pre-tensioning length Al 0 as a constant or posture-dependent parameter, 
and assumed that only the muscle stiffness is under direct neural control. In future extensions of this work, 
we will explore the case where the pre-tensioning values are also regulated by neural drive and consequently 
change the muscle force. 
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